MagnetoResistance of graphene-based spin valves. 
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We study the magnetoresistance of spin-valve devices using graphene as a non-magnetic material 
to connect ferromagnetic leads. As a preliminary step we first study the conductivity of a graphene 
strip connected to metallic contacts for a variety of lead parameters, and demonstrate that the 
resulting conductivity is rather insensitive to them. We then compute the conductivity of the 
spin-valve device in the parallel and antiparallel spin polarization configurations, and find that it 
depends only weakly on the relative spin orientations of the leads, so that the magnetoresistance 
MR of the system is very small. The smallness of MR is a consequence of the near independence 
of the graphene conductivity from the electronic details of the leads. Our results indicate that, 
although graphene has properties that make it attractive for spintronic devices, the performance of 
an graphene-based spin-valve is likely to be poor. 

PACS numbers: 73.20.-r,73.22.-f,72.25.-b 
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I. INTRODUCTION. 

Recently it has become possible to isolate an individual 
graphene layer [l[, a two dimensional crystal of carbon 
atoms packed in a honeycomb lattice. When deposited on 
top of a doped dielectric substrate, the density of carriers 
moving in the graphene sheet can be modified by apply- 
ing a gate voltage [H, At low energies, carriers moving 
in a graphene sheet obey the Dirac equation, so that 
graphene offers the interesting possibility of studying the 
properties of Dirac fermions. Apart from the interest- 
ing fundamental physics of this new system, graphene 
is attracting attention as a promising new material for 
microelectronic applications. 

The conductivity of graphene tends to a minimum 
value when the density of extra carriers tends to zero 
[3, S3- Theories predict a universal value, a u = 4e 2 /nh, 
for the conductivity of a homogeneous and impurity-free 
undoped graphene sheet [E IS H H, II] • However experi- 
mental values of the minimum value of the conductivity 
in graphene sheets is between five and ten times larger 
than the theoretical prediction. This discrepancy may 
occur because in neutral graphene the system breaks up 
into electron and hole puddles [l(|, so that transport 
could occur through the resulting hole-like and electron- 
like regions [ll[ . The charge puddles may appear in 
order to screen impurities which are invariably present 
12 1 , or may be induced by ripples in the graphene sheets 

ii q m m, E3. 

In the ballistic regime, theoretical work fiH [T^. I20I l2l| 
has shown that the conductivity of an intrinsic graphene 
sheet of width W and length L takes the universal value 
in the W/L ~ * oo limit. Such a ballistic approximation 
is valid when the mean free path of the carriers is larger 
than the sample dimensions. This was the case, for exam- 
ple, for devices used in Ref. l22|, which confirmed that the 
conductivity of wide and short graphene ribbons tends to 
the universal value 4e 2 /irh. 



Graphene exhibits room temperature mobilities above 
10 5 cm 2 jVs 0, implying that electrons in graphene 
sheets can move very long distances without scatter- 
ing. For short range scatters the mean free path can 
be as large as lOOOnm (23[. From these results we expect 
the ballistic approximation to be ap pro priate to describe 
transport in graphene nanoribbons [24j, |25| | . In addition 
the small spin-orbit coupling of carbon atoms [26] results 
in a long spin lifetime for carriers in graphene sheets. 
This makes graphene a very good candidate for micro- 
electronic and spintronic applications. 

In developing new spintronic devices, it is very im- 
portant to find non-magnetic materials where a spin- 
polarized current can be injected and flow without be- 
coming depolarized. The most popular existing spin- 
tronic devices are spin valves. These devices use the fact 
[27l | that the electrical resistance of a non-magnetic ma- 
terial connected to spin-polarized source and drain leads 
depends strongly on their relative spin-orientation. Spin 
valves are promising candidates for systems that may 
transform spin information into electrical signals. These 
devices perform best when the spin relaxation time of 
the non-magnetic material is long, making graphene a 
good candidate for this component of a spin- valve. More- 
over, the combination of weak spin-orbit coupling and low 
hyperfine interaction of the electron spin with the car- 
bon nucleus makes graphene a good candidate for other 
spintronic applications such as spin qubits [Hj], three- 
terminal devices [IE H3] and spin filters 31 1 . 

Recently several groups have performed non-local four- 
probe measurements [H, HE H3 of graphene, connected 
to ferromagnetic electrodes, and have demonstrated the 
presence of spin currents between injector and detector. 
These experiments provide proof in principle of the pos- 
sibility of injecting a spin current into graphene, with a 
spin relaxation length larger than 2/im at room tempera- 
ture (32^ . Moreover, a magnetoresistance of several hun- 
dred Ohms was observed in a spin valve where graphene is 
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FIG. 1: (Color online) Schematic representation of the spin 
valves studied in this work. We plot the density of states of 
the central graphene region and of the ferromagnetic source 
and drain leads. In (a) the leads are single orbital band met- 
als where the center of the minority spin band is shifted with 
respect to the majority spin band, (b) corresponds to the 
case of a ferromagnetic transition metal with a spin polarized 
narrow d band and a paramagnetic wide s band. The hori- 
zontal dashed line indicates the chemical potential jj,. In this 
schematic picture y, — 0. In both cases, the spin polarizations 
of the leads is illustrated in the parallel configuration. For 
the antiparallel configuration the spin labels of each band in 
the drain lead should be reversed. 

contacted by two ferromagnetic electrodes 35] . These ex- 
periments clearly demonstrate the potential for graphene 
in spintronic devices. 

In this work we study the magnetoresistance of a 
graphene-based spin valve. This is a three component 
device, with a first ferromagnetic lead used as a spin po- 
larizer, a non-magnetic spacer - graphene in our case 
- and a second ferromagnetic lead used as analyzer. We 
consider two kinds of electrodes (see Fig. [TJ. (a) A single 
orbital band metal where the center of the minority spin 
band is shifted with respect to the majority spin band in 
such a way that the material is ferromagnetic. This is 
a simplified version of a dilute magnetic semiconductor, 
(b) A ferromagnetic transition metal (e.g., cobalt) with 
a spin-polarized narrow d band and a wide paramagnetic 
s band at the Fermi energy. 

We assume the carrier mean free path in graphene is 
longer than the dimension of the graphene part of the de- 
vice, so that the transport properties can be calculated in 
the ballistic approximation. In addition, since the spin- 
orbit coupling in graphene is very small, it is appropriate 
to assume that the spin diffusion length is sufficiently 
long that the carriers do not undergo spin flips while 
traversing the graphene. Therefore we model the trans- 
port in terms of two independent spin channels [36| . In 
this work we will call this the independent current model. 

This paper is organized as follows. In Section II we 
describe the Green's function formalism for computing 



the conductance of a ballistic system [37} , and show how 
this leads to a conductivity for undoped graphene. In 
Section III we then apply the formalism to compute the 
magnetoresistance of a wide piece of graphene contacted 
by two single orbital band ferromagnetic leads. In Section 
IV we discuss the case of transition metal leads, modeled 
as conductors with two orbital bands, one narrow and 
one wide. We then compute the magnetoresistance of 
this system. Finally, we summarize our results in Section 
V. 



II. CONDUCTANCE AND CONDUCTIVITY OF 
A GRAPHENE STRIP 

We begin by first calculating the conductance of a 
graphene strip connected to single orbital metallic leads 
[see Fig. (Ha)] for an arbitrary alignment of the lead band 
bottoms and the Dirac points of the graphene system. 



A. System description. 

We consider a stripe of graphene as illustrated in Fig. 
(5][a), with lattice parameter a, attached to metallic leads 
which are each modelled by a commensurably matched 
square lattice with the same lattice parameter [2(| [2l[ . 
The graphene strip has length L in the x direction and 
is infinitely wide in the y direction, i.e., W — > oo. With 
this geometry, we can define a unit cell, infinitely long 
in the x direction, that is periodically repeated in the 
y-direction [Fig. H[b)]. In order to match the graphene 
sample to the square lattice leads the number of carbon 
atoms in the unit cell, N, should be a multiple of four, so 
that the length of the graphene sample is L = N/4y/3a. 



B. Hamiltonian 

The electronic properties of the graphene are described 
by a tight-binding Hamiltonian with nearest neighbor 
hopping t and zero on-site energy. The left (L) and right 
(R) metallic leads are also described by tight-binding 
Hamiltonians, with nearest neighbor hopping amplitudes 
tL and tji and onsite energies El and er, respectively. 
The graphene sample is connected to the left and right 
leads with hopping amplitudes t' L and t' R respectively. 

Because of the discrete translational invariance in the 
y direction, we can label the eigenvalues and eigenvectors 
of the Hamiltonian by a momentum k y in the y direction 
that takes values in the range —n/a < k y < n/a. The 
limit W — ► oo is obtained by treating this as a contin- 
uous variable. In order to simplify notation, in the rest 
of the paper the wavevector k y will be given in units of 
l/a. With this the electronic properties of the two di- 
mensional problem is reduced to a set of fc^-dependent 
one dimensional Hamiltonians with on site energies (see 
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FIG. 2: (Color online) (a) Structure of the spin valve stud- 
ied in this work. A graphene sample is connected to metallic 
leads which are modeled as commensurably matched square 
lattices, (b) Unit cell that it is periodically repeated in the y 
direction in the case of a graphene strip connected to sin- 
gle orbital metallic leads, (c) Unit cell in the case of a 
graphene strip connected to ferromagnetic transition metal 
leads. Square (blue) points describe the s atoms and circle 
(red) points represent d atoms. 
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ej = £l + 2tL cos ky 
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and < i < N - 1 
and < i < N - 1 

(2) 



The conductance per spin channel of the system takes 
the form 



(3) 



where /i is the chemical potential and g 1D (E, e ) is the 
conductance of the one dimensional chain with an on- 
site energy eo in the left electrode, evaluated for Fermi 
energy E. (Note that g 1D implicitly depends on the on- 
site energy in the right lead as well.) The one dimensional 
conductance is given by [33] 



g (jjt, 2t L cos k y + e L ) 



T^T^G^Ml 2 (4) 



where T^E) = i(E x - and T N (E) = i(£ N - E^) 
are the injection ratios, and Ei and E^v are the self- 
energy terms due to the leads attached at sites 1 and N 
respectively [37l |. 
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In EqJ3|(7i 1 jv(A0 ^ s ^be effective Green's function for the 
system, evaluated in the graphene region: 



/x - H c - Ei(/i) - Eat(^) ' 



(7) 



Here H c is the one dimensional fc. y -dependent Hamilto- 
nian of the graphene sample described by the on-site en- 
ergies and hopping amplitudes given in Eqs. [T]and[2l 

Finally, the conductivity of the device is related to the 
conductance through geometrical factors, 



W 



(8) 



In two dimensions, conductivity and conductance have 
the same units. However, the former is only useful in 
systems where its value is sample-size independent for 
large L and W. Usually this is only the case for diffusive 
systems. Remarkably, as has been pointed out previously 
[l8L [20} and as we demonstrate explicitly below, for un- 
doped graphene it is true in the ballistic limit, provided 
we take W — > 00 before taking L — ► 00. 

In Fig. [3] we plot G(fi)xL/W in units of the universal 
conductivity for two sets of parameters for the leads, in 
the limit of large W, and for two different values of L. 
One may observe from these results that except for (i = 
the conductance is essentially independent of L. Inter- 
estingly, the conductance is relatively insensitive to the 
details of the metallic leads; its overall shape is mostly 
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FIG. 3: (Color online) Conductivity in units of <r u = 
4/n(e 2 /h), as a function of the chemical potential for two 
different graphene sample lengths and two sets of parameters 
describing the lead band structures. The conductivity shows 
a minimum at the Dirac point which is independent of the 
sample length. 



determined by the graphene density of states. In particu- 
lar the roughly linear rise as fj, moves away from zero may 
be understood as reflecting the linear density of states of 
graphene. The shoulders in the conductance riding on 
this linear background correspond to Fabry-Perot reso- 
nances due to the finite length of the graphene sample 
[38| . The amplitude of these resonances is determined by 
the parameters of the leads. 

In the case of intrinsic (i.e., undoped) graphene - ^=0 
- the electronic transport behaves as if it is diffusive, 
and the conductivity is well-defined. For tj, = tn = t' L = 
t' R = t and el = £r : = we obtain a = \/3/2(e 2 /h), 
in agreement with Ref. I20I . This value changes when 
the centers of the energy bands in the leads change. The 
conductivity reaches the universal value, a u = 4/ir(e 2 /h), 
when el — and er — £r [H, H3|. 

The fi = case may be worked essentially analytically, 
because transport is dominated by values of k y in the 
vicinity of the Dirac points, (0,±|^). In this case and 
for large values of L it is possible to explicitly evaluate 
the transmission probability associated with Hq for each 
valley and spin channel using a transfer matrix method, 
with the result 



4£t 2 sin fci sin kn 
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In this expression the parameter £ — ' ' 1 1 1 
tains the information about the coupling to the leads, the 
wavevector k y is defined with respect the Dirac point, 
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FIG. 4: (Color online) Conductivity evaluated at the Dirac 
point, fi — 0, in units of o u = 4/n(e 2 /h), as a function of the 
source and drain leads parameters. 



and k^ and k R are defined in Eq. [S] For \i at the Dirac 
points there are no individual propagating modes in the 
graphene sample connecting the metallic leads at any fi- 
nite value of k y , and the transmission probability van- 
ishes in the limit L — > oo. 

By integrating the transmission amplitude (Eqj9]) over 
k y we obtain the fi — conductance, 



G- 



^9s9v W sin kh sin kn 



arctan 



/ sin (k L +k B ) \ 
h 7T L sin (fcL+fcij) " " \ I— cos (fci+fcij)/ 

(10) 

where g s =2 and g v =2 are the spin and valley degenera- 
cies respectively. Remarkably, although all the modes are 
evanescent, the total conductance falls off only as l/L, 
so that the electrical transport is Ohmic. This behavior 
arises because the effective length scale of the evanes- 
cent modes are each l/k v , which is arbitrarily large as 
k y — > 0, resulting in the relatively weak l/L behavior 
for G. (Alternatively, the finite value of the transmission 
probability near the Dirac points may be understood in 
terms of virtual electron-hole pair excitations near zero 
energy in the graphene region [181].) By contrast, for sys- 
tems in which /i may lie in a gap, the effective length scale 
for wavefunctions is bounded by a distance that is deter- 
mined by the difference between fj, and the bottom of the 
conduction band for the "conducting" region. Because 
of this maximum length scale G falls off exponentially 
with L, and the conductivity as well as the conductance 
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vanish as L — > oo. Thus clean graphene is rather unique 
in displaying Ohmic behavior. 

Another remarkable feature of this result is that the 
conductivity does not depend on the couplings t' L and 
t' R of the metallic leads with the graphene sample. This 
result is only possible in the limit W — > oo, since for 
any finite value of W the two leads become fully dis- 
connected if either of these parameters vanishes, and the 
conductance must vanish. However, in the infinite width 
limit where the momentum sum becomes an integral, the 
parameter £ through which t' L and t' R enter scales out, 
and the final result is independent of these parameters. 

Finally, it is interesting to examine the dependence of 
the conductivity on the parameters specifying the leads. 
From the definition of fc^ and ka one can see that the con- 
ductance depends on the lead parameters only through 
the combinations (eR—tn)/tR and (e^— In Figure 

H]we plot the conductivity, a = G x L/W, as a function 
of these combinations. The conductivity is maximized 
when tL = Sl and tji — er, independent of the val- 
ues of and tn. Because of conservation of energy and 
of the transverse momentum, k y , transport through the 
graphene ribbon when \i = is possible only when the 
band center of the left and right leads are in the intervals 
2£l < th\ and2iL < th\ respectively. In Ref. I20I 
Schomerus suggests that the conductivity of a graphene 
sample connected to metallic leads is maximum when 
the self-energies of electrons in the leads have the same 
value as the self-energy of bulk graphene with energy 
at the Dirac point (—it). In the square lattice the self- 
energies of the incoming and outcoming electrons, eval- 
uated at the momentum and energy of the Dirac point, 
are —tLe lkL and —tRe lkR , respectively. Our results in- 
dicate that when the real part of the lead self-energies 
= 0) and £_/v(/i — 0) (Eq. O vanish, the conductiv- 
ity of the system is maximized. This condition is similar 
to but less restrictive than that proposed in Ref. [2(J and 
is related to the independence of the graphene universal 
conductivity from the hopping amplitude t. 



III. MAGNETORESISTANCE: SINGLE BAND 
LEADS 

In this section we apply the results derived above to 
find the magnetoresistance of a spin valve device with a 
graphene strip at its center. In a single band ferromag- 
netic metal, the centers of the spin up (eo,t) anc l spin 
down (eo,j.) bands are shifted, as illustrated in Fig. [IJa). 
This implies a relative spin polarization P of the carriers 
at the Fermi energy \x given by 



where K is the complete elliptic integral of the first kind, 
and if and are the hopping matrix elements in the 



spin up and spin down channels respectively. Eq. [TT] is 
obtained directly from the density of states at energy E 
of a square lattice with hopping parameter t and on-site 
energy Eq, which is given by [391 ] 



p(E) = 4t©(2M£-£o|)K 
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The transport through a non-magnetic material con- 
nected to ferromagnetic metals is expected to depend 
strongly on the magnitude and relative orientation of the 
polarizations of the leads. The magnetoresistance (MR) 
is defined as the relative change of the resistance when 
the relative spin orientation of the leads changes from 
parallel, R par a, to antiparallel, R an u- In the ballistic ap- 
proximation the MR can be written as 



MR = 



R„ 



R 



■para 



G 



para 



G„ 



Pa 



G,, 



(13) 



where G para and G an ti are the conductances of the sys- 
tem when the relative polarization of the leads is parallel 
or antiparallel respectively. 

When the non-magnetic material is an insulator the 
transport is through tunneling processes. Assuming that 
tunneling transport is proportional to the the density of 
states at the Fermi energy, Julliere [40j proposed the fol- 
lowing expression for the tunneling MR: 



TMR 



2P 2 



1 + P 2 ' 



(14) 



Jullicrc's expression works rather well for tunneling MR 
even for materials with complicated band structures [4lj ] . 
It predicts correctly that the MR increases as the leads 
become more spin polarized. In this section we study the 
MR when the non-magnetic material is graphene, and 
analyze the results as a function of the spin polarization 
of the leads. 

In the independent current model, and using the bal- 
listic approximation, G par a is the sum of the spin up 
conductance, Eq. [3] evaluated with el = er = £o,T, 
and the spin down conductance, Eq. [3] evaluated with 
El = En — £ o,i- I n the case of antiparallel spin polar- 
ization of the leads G an ti is the sum of Eq. [3] evaluated 
with El = £o,T an( i £ R = e o,i an d Eq- El evaluated with 
e l = £0,1 and e r = e ,T- 

In Fig. [5] we plot the magnetoresistance MR of a 
graphene based spin valve as function of the positions 
of the center of the spin up and spin down bands. We 
consider only intrinsic graphene (/i = 0), and the hop- 
ping parameter in the metallic leads is taken equal to 
that in the graphene part of the device, tL—tR—t. The 
graphene slab is extrapolated to infinite length, although 
for L larger than about 20a the results are essentially the 
same. 

We observe that in general the magnetoresistance is 
small. Only when the parameters are such that one of 
the spin bands is close to or in the forbidden transport 
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FIG. 5: (Color online) (a) Magnetoresistance (MR) and (b) 
spin polarization at the Fermi energy, as a function of the 
centers of the spin bands, of a graphene spin valve. The de- 
vice consists of a long graphene slab separating single band 
ferromagnetic metals, Fig[lja). The chemical potential /^=0, 
and ti, =tR = t. 



region, 2t < |e U(T — t\ does one find MR to be a significant 
fraction of one, and when one of the spin bands is in the 
forbidden region then MR of course reaches its maximum 
possible value. Thus in order to get a moderate value of 
the magnetoresistance a large shift between the center of 
the spin bands in the leads is needed. 

The reason for the smallness of M R is the weak depen- 
dence of the graphene conductance on the parameters of 
the leads [see Figs. [3] and SJ] This weak dependence, 
particularly on the density of states of the incoming elec- 
trons, also implies an absence of any strong correlation 
between the spin polarization of the incoming electrons 
[Fig. HJb)] and the magnetoresistance [Fig. |5ja)]. By 
comparing the values of the polarization and the magne- 
toresistance we observe that the Julliere expression, Eq. 
RBI fails to describe the magnetotransport properties of 
graphene-based transistors. 

Note also that because the conductivity of undoped 
graphene is independent of the tunneling amplitude con- 
necting the graphene with the ferromagnetic metallic 
leads (Eq. fTU|) . the magnetoresistance is also unaffected 
by changes in the quality of the connection between the 
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FIG. 6: (Color online) Spin polarization, evaluated at the 
Fermi energy, as a function of the position along the graphene 
slab of a graphene-based spin valve. We plot the u — case, 
(a) and (b) correspond to antiparallel and parallel configura- 
tions respectively. Square and dot symbols indicate the spin 
polarization on the two types of lattice sites of the honey- 
comb lattice. The band structure parameters are indicated 
in the inset of the figure. The length of the graphene slab is 
L = 20 x V3o. 



leads and the graphene. This indicates that the smallness 
of the magnetoresistance in graphene-based devices is not 
due to a conductivity mismatch, as is the case in metal- 
semiconductor diffusive junctions [42l |. but rather is due 
to the universal minimum conductivity of graphene, so 
that a "shutoff" of one spin channel is difficult to achieve. 

In Fig. [6] we have plotted the spin polarization as a 
function of position in our graphene-based spin valve. 
For this calculation the lead parameters were chosen so 
that the spin polarization in the leads was approximately 
60%. Several points are worth noting. (1) The spin po- 
larization decays to a zero over a length scale of about 20 
lattice parameters. (2) The spin polarization is oriented 
in opposite directions for electrons on different sublat- 
tices. (3) The total induced spin in the graphene region 
vanishes. These effects are explained by the strong ten- 
dency for local magnetic moments in graphene to orient 
ferromagnetically for sites on the same sublattice, and an- 
tiferromagnetically for sites on different sublattices [43j j - 
We note that this result does depend on the fact that in 
our chosen geometry, zigzag edges present themselves to 
the leads, so that all the graphene sites contacting the 
leads are on the same sublattice. The sensitivity of the 
electronic states to the edge geometry is a well-known 
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FIG. 7: (Color online) (a) Magnetoresistance (MR) and 
(b)spin polarization at the Fermi energy as function of 
the chemical potential in the graphene layer and the spin 
down band center. The center of the spin up band is lo- 
cated in the optimal conductivity position, so,i=t. We take 
tL~tR—t' L —t' R . The length of the graphene slab is set to 
L = 20 x \^3a. For larger values of L the results are essen- 
tially identical. 



IV. MAGNETORESISTANCE: 
FERROMAGNETIC TRANSISTION METAL 
LEADS 

Energy bands for ferromagnetic transition metals can 
be described approximately by using two bands, a d-band 
characterized by a width 4t<j and center position Ed, and 
an s-band with width At s and center at e s (47[. Because 
the shape of the atomic orbitals, the s-band is much wider 
that the d-band. Cobalt, for example, is a transition 
metal with a minority spin d-band shifted up in energy 
relative to a majority spin d-band. The s-band is nearly 
spin-unpolarized, and the Fermi energy crosses both the 
majority spin d-band and the unpolarized s-band, as il- 
lustrated in FigJTJb). For studying the transport proper- 
ties of such systems, it is a good approximation to neglect 
the minority spin d-band located at high energy. 

We wish to consider the same device as studied in Sec- 
tion III with leads characterized by bandstructures of this 
sort. This can be modeled by assigning two orbitals to 
each site in the leads, one for the s-band and the other 
for the d-band. We can then assign on-site energies e Sj d 
respectively for the two orbitals, and hopping matrix el- 
ements t s and td- These orbitals then connect to the p z 
carbon graphene orbitals through hopping amplitudes t' s 
and t' d respectively. Although cobalt is a three dimen- 
sional metal, for simplicity we treat it here as two di- 
mensional. The resulting one-dimensional problem that 
one needs to treat for each k y is illustrated in Fig. [UJc). 
Finally, for this section we will restrict our discussion to 
the case of intrinsic graphene, /i = 0. 



Conductivities 



property of graphene 0, |H, HgJ 

Finally, we analyze the magnetoresistance for doped 
graphene in the three stripe geometry. In this case the 
transport through the graphene part of the device is bal- 
listic in the usual sense, and the MR has to be calculated 
as the relative change in the conductance, Eq. [T3J In Fig. 
[TJa) we plot the MR for this case, as a function of the 
chemical potential in the graphene, and of the center of 
the spin down band. In addition we plot the spin polar- 
ization at the leads. As in the case of undoped graphene 
the MR is small (although not quite so small as in the 
undoped case), and only moderate values of MR result 
from large shifts between the centers of the spin up and 
spin down bands. From comparing the two panels of Fig. 
El we again observe that the values of MR and the spin 
polarization at the leads are largely uncorrelated. 



As there are two band orbitals in the transition metal 
leads, in order to compute the conductivity with different 
spin polarization orientations in the leads we will need 
to evaluate several partial conductivities, representing 
transmission between different combinations of orbitals. 
These are: 

• The conductivity across a graphene strip attached 
to single band spin-polarized metals, with the 
bands in the two leads centered at the same en- 
ergies and having equal bandwidths, 

<r v , v = ^--G(k u ), (15) 

where G(k u ) is the conductance (Eq. [TO)) evalu- 
ated with kr J —kfi=k L , — arccos( t "^ e " )■ This con- 
tribution enters when we consider the transmission 
of spin-minority electrons in the parallel configu- 
ration, so that the d-bands of the minority spin 
are high in energy in both leads and are irrelevant. 
Thus only the band index v = s for this contribu- 
tion will be relevant in our calculation. 
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The conductivity of the device where the source and 
drain metals consist of two conducting bands, s and 
d. In this case the conductivity per spin channel is 



0~sd,sd 



— (t s sin k s + td sin k d ) 2 x 
a 7r 



arctan 



v /4| n |4-(a-H a « 2 )2 
2|a| 2 — a 2 — a* 2 ' 



V4|a| 4 - (a 2 + a* 2 ) 2 



with 



(16) 



(17) 



This contribution enters for spin-majority electrons 
when the lead polarizations are parallel. 

• Finally it is necessary to evaluate the conductivity 
in the case where the only available band in the 
source metal is the s-band, whereas both bands, s 
and d, are available in the drain lead, or vice- versa. 
These two situations arise when the lead polariza- 
tions are antiparallel. The required conductivities 
are both 



?s,sd — 0~sd : s 



e ^9v ( 2 ■ 2 j 

= t [t s an k s 



t s t d sin k d sin k s 



arctan 



y/Hb^-i^+b' 2 )- 
2\b\ 2 -b ,2 -b* 2 



with 



V4|6| 4 - (b 2 + b* 2 ) 2 



(i s e ° +t d e lkd )e lk °t 



(18) 



(19) 



The above conductivities are independent of the values 
of the contact hopping amplitudes between the graphene 
orbitals and the s and d orbitals of the transition metal 
leads. 

A very interesting result which emerges from Eqs. 
(fTS)) . (fTo]) and (fT5)) is that the conductivity of intrinsic 
graphene is nearly independent of the number of bands in 
the drain and source leads. For example, if we take e s =t s 
and Sd—td we find that all the conductivities are the 

same: cr ss =add=Vsd,s=Vsd,sd = (Note that we 

are evaluating the conductivity per spin channel). Simi- 
larly, if we center all the bands at zero energy, e s =e d =0, 
all the conductivities take the value V3/4^-, indepen- 
dent of the hopping parameters. Fig. [8] illustrates this 
behavior over a range of lead parameters. These results 
reflect the remarkably weak dependence of the graphene 
conductivity on the electronic structure of the contacts, 
and in particular on the density of states of the metallic 
leads at the Fermi energy. 



B. Magnetoresistance 

For computing the magnetoresistance we have to eval- 
uate the conductivity in the parallel and antiparallel spin 
polarization configurations of the leads. 




| 0.92 
O 0.90 



0.88 4 



FIG. 8: (Color online) Conductivity per spin channel of a 
a graphene slab connected to metallic leads with two (s and 
d) bands each, o- s d, s d and to a two orbital (s and d) source 
lead and a single orbital (s) drain lead, o s d >s - We plot the 
conductivity as a function of the band center of the d-band. 
We assume that the s-band is much wider than the d band, 
taking t„=5t and td=0-2t. Plotted are the e v =t v and the e„=0 
cases. Note that in all the cases the conductivity depends 
weakly on the parameters of the metallic leads. 



In the parallel configuration, majority spin electrons 
in the s and d bands of the metallic source are injected 
in the graphene slab and received in the majority spin 
s and d bands of the metallic drain. The minority spin 
electrons can go just from the source s-band to the drain 
s-band. Therefore the conductivity is 



0~para — 0~sd,sd i 0~ ss . 



(20) 



In the antiparallel configuration the majority spin elec- 
trons in the source can be in the s or d band. Upon 
passing through to the drain lead these are now minor- 
ity electrons, which can only reside in the s-band. The 
inverse situation occurs for the minority spin carriers in 
the source lead. Thus the conductivity in the antiparallel 
configuration takes the form 



0~anti 



0~sd,. 



(21) 



Since by symmetry o- s-s d~o' s d,s, the magnetoresistance 
takes the form 



MR = 



0~sd,sd + 0- SiS — 2(T S 
0~sd ,sd + °~s ,s 



(22) 



In Figure 02a) we plot our calculated magnetoresis- 
tance for this system. We assume the s-band to be much 
wider than the G?-band, and plot the results as functions 
of the energy center of the c?-band, and for different posi- 
tions of the center of the s-band. In Figure (Hb) we plot 
the corresponding spin polarization of the ferromagnetic 
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leads. From Fig. [9] we immediately see that the magne- 
toresistance is generically very small, which we can again 
understand as a consequence of the relative insensitivity 
of the conductivity to the details of the leads. Inter- 
estingly, this same insensitivity means that there is no 
simple relationship between the value of the spin polar- 
ization in the leads and MR, as can be seen by comparing 
Figs. [HJa) and[5]Jb). (Note that the peak appearing in 
the polarization at Ed = is due to a van Hove singu- 
larity in density of states for a two dimensional square 
lattice in the tight binding Hamiltonian.) 

Finally, we note that recent experiments [35] on a 
graphene-based spin valve did lead to a significant MR 
for undoped graphene. Our results show that this result 
cannot be understood purely on the basis of a clean, non- 
interacting electron model. However, the ferromagnetic 
contacts reported by Hill et al. in Ref. HH exhibited 
strongly non-Ohmic, nonlinear I-V behavior which one 
would not usually associate with normal metallic con- 
tacts 48]. When similar devices with Ohmic contacts 
(using a Ti sublayer) were fabricated, the magnetoresis- 
tance previously observed completely disappeared [48| . 
These experimental findings support our results that spin 
valves employing clean graphene and simple, homoge- 
neously magnetized leads should have relatively very low 
MR. 



V. SUMMARY AND CONCLUSION 

In this work we undertook a detailed study of the con- 
duction properties of wide graphene strips, with two dif- 



ferent models for the source and drain leads. We recon- 
firmed that for undoped graphene, the system can be de- 
scribed by a conductivity in the L — > oo limit even when 
defects are absent from the system, and examined this 
behavior with respect to a broad range of lead parame- 
ters. The resulting conductance turns out to be relatively 
insensitive to these. 



We then used these results to compute the magnetore- 
sistance of a simple three stripe spin-valve device with 
graphene acting as the non-magnetic material between 
the ferromagnetic leads. Two types of ferromagnetic lead 
systems were considered: one with a single (s) orbital for 
each spin state, with band centers separated in energy 
to induce spin polarization, and another with a narrow 
d band which was taken to be spin-polarized. It was 
found that the magnetoresistance MR was rather small 
for most circumstances in both cases, largely due to the 
insensitivity of the conductivity with respect to condi- 
tions in the leads. 
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FIG. 9: (Color online) (a) Magnetoresistance of a graphene 
strip connected to ferromagnetic transition metal leads, (b) 
Spin polarization at the Fermi energy of the ferromagnetic 
leads. The d band is assumed to be much wider than the s 
band, t s =5t and td=0.2t, with t the graphene hopping pa- 
rameter. We present the results for different positions of the 
center of the s band. 
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